DD Correlations as a Sensitive Probe for Thermalization in High Energy Nuclear Collisions 
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We propose to measure azimuthal correlations of heavy-flavor hadrons to address the status of thermalization 
at the partonic stage of light quarks and gluons in high-energy nuclear collisions. In particular, we show that 
hadronic interactions at the late stage cannot significantly disturb the initial back-to-back azimuthal correlations 
of DD pairs. Thus, a decrease or the complete absence of these initial correlations does indicate frequent 
interactions of heavy-flavor quarks and also light partons in the partonic stage, which are essential for the early 
thermalization of light partons. 

PACS numbers: 25.75.-q 



Introduction 

Lattice QCD calculations, at vanishing or finite net-baryon 
density, predict a cross-over transition from the decon- 
fined thermalized partonic matter (the Quark Gluon Plasma, 
QGP) to hadronic matter at a critical temperature Tc ~ 150- 
180MeV[l]. 

Measurements of hadron yields in the intermediate and high 
transverse momentum (pj) region indicate that dense matter 
is produced in AuH-Au collisions at RHIC ||2|, Ij]. The ex- 
perimentally observed large amount of elliptic flow of multi- 
strange hadrons l^, such as the (|) meson and the Q. baryon, 
suggest that collective motion develops in the early partonic 
stage of the matter produced in these collisions. A crucial is- 
sue to be addressed next is the thermalization status of this 
partonic matter. 

Heavy-flavor (c, b) quarks are particularly excellent 
tools [5, 6, 7] to study the thermalization of the initially cre- 
ated matter. As shown in Fig.[Tl their large masses are almost 
exclusively generated through their coupling to the Higgs field 
in the electro-weak sector, while masses of light quarks (u, 
d, s) are dominated by spontaneous breaking of chiral sym- 
metry in QCD. This means that in a QGP, where chiral sym- 
metry might be restored, light quarks are left with their bare 
current masses while heavy-flavor quarks remain heavy. Due 
to their large masses (^ ^qcd), the heavy quarks can only 
be pair-created in early stage pQCD processes. Furthermore, 
their production cross sections in nuclear collisions are found 
to scale with the number of binary nucleon-nucleon colli- 
sions ISl |9|]. In the subsequent evolution of the medium, 
the number of heavy quarks is conserved because the typical 
temperature of the medium is much smaller than the heavy 
quark (c, b) masses, resulting in negligible secondary pair 
production. In addition, the heavy quarks live much longer 
than the lifetime of the formed high-density medium, de- 
caying well outside. These heavy quarks (c, b) can partic- 
ipate in collective motion or even kinetic ally equilibrate if, 
and only if, interactions at the partonic level occur at high 



frequency. The idea of statistical hadronization of kineti- 
cally equilibrated charm quarks [10] even predicted signifi- 
cant changes in hidden charm hadron production 1 1 11 . Hence, 
heavy-flavor quarks are an ideal probe to study early dynamics 
in high-energy nuclear collisions. 
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FIG. 1: (Color online) Quark masses in the QCD vacuum and the 
Higgs vacuum. A large fraction of the light quark masses is due to 
the chiral symmetry breakin g in the QCD vacuum. The numerical 
values were taken from Ref. 1131. 



Recent STAR and PHENIX results on elliptic flow and 
nuclear modification factors of non-photonic electrons indi- 
cate that charm quarks might indeed participate in the col- 
lective motion of the matter produced in AuH-Au collisions at 
RHIC. To explain the data, a large drag diffusion coefficient 
in the Langevin equation, describing the propagation of charm 
quarks in the medium, is found to be necessary 11311 . Two- and 
three -body interactions lll4lll5ll of heavy-quarks and resonant 



rescattering lllTll in the partonic stage seem to become impor- 
tant. These investigations suggest that heavy-quarks actively 
participate in the partonic stage. 

In this paper, we study the change of azimuthal correla- 
tions of D and D meson pairs as a sensitive indicator of fre- 
quent occurrences of partonic scattering. Since heavy-flavor 
quarks are pair-created by initial scattering processes, such as 
gg — > cc, each quark-antiquarkpair is correlated in relative az- 
imuth A(|) due to momentum conservation. In elementary col- 
lisions, these correlations survive the fragmentation process to 
a large extent and hence are observable in the distribution of 
the relative azimuth of pairs of D and D mesons. We evalu- 
ate by how much these correlations should be affected by the 
early QGP stage and by the hadronic scattering processes in 
the late hadronic stage. 

Results and Discussions 

We start by reviewing the D and D angular correlations in 
pp collisions. The Monte Carlo event generator PYTHIA 1 18] 
reproduces well the experimentally observed correlations of 
D mesons, measured at fixed target energies 111911 . Figure |2] 
shows such correlations as calculated by PYTHIA (v. 6.139) 
for pp colhsions at -/s = 200 GeV. The A(|) distribution is 
peaked at 1 80°, especially for high pj D mesons, as one would 
expect for back-to-back pairs stemming from hard scatterings 
of partons. 
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FIG. 2: (Color online) Correlations in relative azimuth, A(|), of DD 
pairs from pp collisions at yG = 200 GeV, as calculated by PYTHIA 
(v. 6.139) for different ranges of single D meson pj. 



To explore how the QCD medium generated in central ultra- 
relativistic nucleus-nucleus collisions influences the correla- 
tions of D and D meson pairs, we employ a non-relativistic 
Langevin approach which describes the random walk of 
charm quarks in a QGP and was first described in Refs. [5, 



^ = -y{T)p + l\, 



(1) 



where fj is a Gaussian noise variable, normalized such that 
{r[i{t)^j{t')) — a{T)5ij5{t — t') with i,j indexing directions. 
Both the drag coefficient y and the momentum-space diffu- 
sion coefficient a depend on the local temperature, T. We 
use the same parameterization for y as in Ref. ||5|], y(r) = 
aT^, with fl = 2 • 10"^ (fm/c)"' MeV"^ and also neglect a 
momentum dependence of y. y and a are related by the 
fluctuation-dissipation relation in equilibrium, {pj) = a/2y. 
As in Ref. [5], a is calculated from the above equation with 
(/??) — 1 .33 nic T, with a charm quark mass jtic = 1 .5 GeV/c^. 

For simplicity, the initial conditions for central collisions 
of heavy nuclei (Au or Pb) are the same as in Ref. ||51]. We 
assume a plasma in thermal equilibrium occupying a cylin- 
der of fixed radius R, equal to the radius of the colliding nu- 
clei (R = 7 fm in the present calculations), at the time of its 
full formation Tq- After that, the plasma evolves according to 
Bjorken's hydrodynamics, with a temperature dropping like 
T = Tq (tq/t)^'^. Before the time Tq, the plasma is changing 
rapidly and is not fully equilibrated. We assume (as done in 
Ref. ||5|]) that the plasma is already in equilibrium before Tq, 
with T = Tq and with the charmed quarks diffusing in the same 
way as they do after Tq. In order to isolate the effects purely 
due to parton-parton re-scattering in the outlined medium, we 
generated c and c quarks with the same pi and zero longi- 
tudinal momentum back-to-back, i.e. with A(|) = 7t at a time 
of the order of I /m^ ~ 0. 1 fm/c, and used a delta function 
for fragmenting the charm quark into a charmed hadron at the 
hadronization stage. The radius r where each pair is created 
is randomly generated from a distribution reflecting the num- 
ber of binary nucleon-nucleon collisions taking place at that 
radius, p{r)dr°^ {R^ — r^) 27trdr. The angular distribution of 
the initially produced charm quarks in the transverse plane is 
isotropic. 

The evolution of the charm momenta with the Langevin 
equation is stopped when T reaches the critical temperature 
Tc = 165 MeV or when the charm quark leaves the QGP vol- 
ume. Furthermore, to get a first estimation of the QGP ef- 
fect on the charm quark pairs azimuthal correlation, we omit 
the possible contribution of the mixed phase. Figure [3](a) 
shows the results for the D meson (charm quark) pairs an- 
gular correlations for different initial charm quark pj values, 
for Tq ^ 300 MeV and To = 0.5 fm/c (typical values for RHIC 
collisions). We see that the fastest charm quarks (represented 
by the pj = 3 GeV/c red curve) are able to escape from the 
QGP without suffering significant medium effects, while the 
slower quarks (see the pj — 0.5 GeV/c black line) have their 
pair azimuthal correlation almost completely smeared out by 
the interactions in the medium. 

Figure [3](b) shows the corresponding results for Tq = 
700 MeV and To — 0.2 fm/c, values representative of LHC 
energies. Here, the interactions of the charm quarks with the 
medium are so frequent that only the most energetic charm 



quarks preserve part of their initial angular correlation; low 
Pj pairs can even be completely stopped by the medium. Be- 
cause the c and the c quarks of a given pair are created to- 
gether, in the same space point, the pair has a higher escape 
probability if both quarks escape the medium from the side 
of the fireball where the thickness is smaller. Thus, the A(|) 
distribution is shifted towards smaller values. Presently, there 
is only longitudinal flow in the hydro-dynamical model used 
in the calculations. We expect that the strong radial flow will 
further enhance the same side correlations in the case of fully 
kinetically thermalized charm pairs. Because the two quarks 
of a pair are created in the same position in the medium, they 
will be pushed in the same direction by the radial flow. In 
fact, this mechanism is a known non-flow effect in the mea- 



surement of collective flow 112111 . 

The above results were obtained with a drag coefficient esti- 
mated with perturbative QCD adopting a large coupling con- 
stant a, = 0.6 1 20]. The recent pQCD calculations 111 [ntl 
show a factor of 2-3 smaller drag coefficient. Besides, as 
shown in Ref. (iT^, non-perturbative contributions that arise 
from quasi-hadronic bound states in the QGP might be impor- 
tant. The presence of these resonances at moderate QGP tem- 
peratures would result in a much larger drag coefficient. Since 
exact values of the drag coefficient from lattice QCD calcula- 
tions do not exist, we show the sensitivity of our calculations 
when varying the parameter a within reasonable values. Fig- 
ures [3](c) and (d) show the drag coefficient dependence of the 
cc angular correlations for high-energy c quarks: pj = 3 and 
10 GeV/c for Tq = 300 and 700 MeV, respectively. The cc 
angular correlations are washed out when a is increased by 
around a factor of 5 in the first case and by around a factor of 
2 in the second, with respect to the pQCD value. 
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FIG. 3: (Color online) Correlations in relative azimuth NS) of DD 
pairs from Langevin calculations with Tq = 300 MeV (RHIC) and 
700 MeV (LHC). The upper part shows the dependence of the corre- 
lations on the initial pi of the c quark; the lower part shows the drag 
coefficient dependence. 



While the initial correlations of cc pairs are clearly affected 
by the formed QCD medium, and sensitive to its temperature 
and drag coefficient, it is important to evaluate to which ex- 
tent the loss of correlations can be mimicked by hadronic in- 
teractions. Therefore, we have investigated by how much the 
cc correlations deteriorate due to hadronic scattering in the 
late stage of high-energy nuclear collisions. For this purpose, 
we used the microscopic hadronic transport approach UrQMD 
V.2.21I2II2I. 

This microscopic transport approach is based on the covari- 
ant propagation of constituent quarks and diquarks accompa- 
nied by mesonic and baryonic degrees of freedom. It sim- 
ulates multiple interactions of in-going and newly produced 
particles, the excitation and fragmentation of color strings, 
and the formation and decay of hadronic resonances. A phase 
transition to a QGP state is not incorporated into the model dy- 
namics. In the latest version of the model, PYTHIA (v. 6.139) 
was integrated to describe the energetic primary elementary 
collisions [24]. Measured particle yields and spectra are well 
reproduced by the approach 12411 . 

In this model, D mesons stem from the very early 
stage high-energy nucleon-nucleon collisions, calculated with 
PYTHIA. For their propagation in the hadronic medium, we 
consider elastic scattering of D mesons with all other hadrons. 
The hadronic scattering cross-section for D mesons is gener- 
ally considered to be small [25] and we take 2 mb in our cal- 
culation. To evaluate the sensitivity of the calculations to this 
parameter, we have repeated them using a ten times higher 
cross section, 20 mb. The final state angular distribution of 
the elastic scattering between D mesons and other hadrons can 
affect the DD correlations. An isotropic angular distribution 
in the scattering process, motivated by the idea of interme- 
diate resonances formation (essentially D+Jt ^^D*), has the 
strongest effect on the charm correlations. For comparison, 
we have also used a forward peaked distribution, simulated 
according to the distribution /7(cos9) <=<: exp(7cos9), where 9 
is the scattering angle in the center of mass system of the de- 
caying resonance. 

The results of the UrQMD calculations are shown in Fig.H] 
for minimum bias AuH-Au reactions at y^sj^= 200 GeV. To 
better illustrate the change in the angular correlations of the 
DD pairs, we show the ratio between the final distribution, 
affected by the evolution with UrQMD, and the initial one, 
directly obtained from PYTHIA (see Fig. |2]i. For this spe- 
cific analysis, all DD pairs were selected, irrespective of the 
p-\ of the D mesons. We observe that the DD pair correlation, 
evaluated through the relative azimuthal angle distribution, is 
barely affected by hadronic interactions. Except in the some- 
what drastic case of using a 20 mb hadronic cross-section and 
even in this case only if we use an isotropic scattering angular 
distribution (Fig.|4](b)) the initial correlations are weakened. 

We also find that most scatterings happen at the earliest 
times, when the hadron density is very high. At RHIC en- 
ergies, the (pure) hadronic stage is presumably shorter than 
assumed in the present calculations. Therefore, the hadronic 
stage should not induce visible changes to the measured DD 



angular correlations. Thus, a change of the DD angular cor- 
relation in heavy-ion collisions, is dominated by frequent 
parton-parton scatterings occurring in the QGP phase. 
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FIG. 4: Ratios between the A(|) distributions of DD pairs produced in 
minimum bias Au+Au reactions at a/s^= 200 GeV, before and after 
hadronic rescattering, using forward (left panel) and isotropic (right 
panel) angular distributions. 

The calculations reported above show that heavy-flavor cor- 
relations provide a sensitive tool to directly access the sta- 
tus of equilibration of the partonic medium. A similar pic- 
ture emerges from the observation and suppression of back-to- 
back correlations of unidentified high pj hadrons (jets). Both 
features originate from the propagation of partons (heavy 
quarks or jets) in the medium. However, D mesons (and B 
mesons) have the advantage that they can still be identified, 
even if they lose a significant fraction of their energy and get 
kinetically equilibrated. 

Compared to the corresponding balance functions |26] of 
electric charge or strangeness, heavy-flavor correlations have 
obvious advantages: Heavy quarks and anti-quarks are pair 
created in the early stage hard processes only and their sub- 
sequent annihilations are negligible. Hence, both numbers of 
heavy quarks and anti-quarks are conserved separately. On 
the other hand, only the net-charge or total strangeness is con- 
served while the number of positive (and negative) charges or 
strange and anti-strange particles can be changed due to new 
pair production and annihilation throughout the whole life- 
time of the fireball. 

Conclusions and Outlook 

In summary, we argue that the observation of broadened an- 
gular correlations of heavy-flavor hadron pairs in high-energy 
heavy-ion collisions would be an indication of thermalization 
at the partonic stage (among light quarks and gluons). We 
have seen that hadronic interactions at a late stage in the colli- 
sion evolution cannot significantly disturb the azimuthal cor- 
relations of DD pairs. Thus, a visible decrease or the com- 
plete absence of such coiTelations, would indicate frequent 
interactions of heavy-flavor quarks and other light partons 
in the partonic stage, therefore imply early thermalization of 
light quarks (which would be the main scattering centers) in 
nucleus-nucleus collisions at RHIC and LHC. 

These measurements require good statistics of events in 
which both D mesons are cleanly reconstructed. A complete 



reconstruction of the D mesons (i.e. of all their decay prod- 
ucts) in full azimuth is essential to preserve the kinematic in- 
formation and to optimize the acceptance for detecting corre- 
lated D meson pairs. Solid experimental measurements in pp 
and light-ion collisions, at the same energy, are crucial for 
detailed studies of changes in these azimuthal correlations, 
and should be performed as a function of pj. Upgrades of 
the STAR and PHENIX detectors with micro-vertex capabil- 
ities II27II and direct open charm reconstruction should make 
these measurements possible. 

Note that more than one pair of DD could be produced in 
high-energy nuclear collisions. At mid-rapidity, the total num- 
ber of charm pairs is in the order of 2-4 and 6- 1 at RHIC and 
LHC, respectively. The multiple production of un-coiTelated 
D mesons will give an additional background to the measured 
angular correlations. Our tests show that the effect can be 
removed by the mixed event method. The mixed-event sub- 
tracted correlation function can be compared with the calcu- 
lated results. In addition, the next to leading order contribu- 
tions, such as gluon splitting, to charm production may be- 
come more important in higher energy collisions. These DD 
pairs are produced without a clear back-to-back correlations 
and will lead to a much weaker peak at 180 degrees. In order 
to understand the thermalization at RHIC, more theoretical 
analyses in this direction are called for. 
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